Admixtures to d-wave gap symmetry in untwinned YBa2Cu307 superconducting films 

measured by angle-resolved electron tunneling 
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We report on an afe-anisotropy of Jc\\b/Jc 



1.8 and IcRn\\b/ IcRn\\a — 1-2 in ramp-edge junc- 



tions between untwinned YBa2Cu307 and s-wave Nb. For these junctions, the angle 9 with the 
YBa2Cu307 crystal 6-axis is varied as a single parameter. The _R„A(&)-dependence presents 2-fold 
symmetry. The minima in IcRn at = 50° suggest a real s-wave subdominant component and 
negligible d^y-^jwe or imaginary s-wave admixtures. The Jc-Rn (6')-dependence is well- fitted by 83% 



15% isotropic s- and 2% anisotropic s-wave order parameter symmetry, consistent with 



At/Aa ?^ 1.5. 



PACS numbers: 74.20.Rp, 74.50.-fr, 74.72. Bk, 74.78.Bz, 85.25.-j 



Phase-sensitive 



experiments 



and tunnel 

spectroscopjiii have provided rich evidence for the 
sign change of the pair wave function in the crystal ab- 
plane of high-T^ superconductors. Insight in the extent 
of subdominant admixtures to the dx2_y2--wave symme- 
try is less vifell established. They are of high importance 
for the basic understanding of high-Tc superconductivity 
and the design of novel d-wave based Josephson devices, 
but also for standard high-Tc junctions. They determine 
for instance the exact position of the nodes, and the 
amount of a5-anisotropy. 

In YBa2Cu307 strong anisotropy in the electronic 
structure has been reported, which can be interpreted 
as an effective mass anisotropy along the a- and b- 
axes: An elongated vortex shape by scanning tunnel- 
ing spectroscopy'' suggests 50% anisotropy. Sixty per- 
cent anisotropy is found in the London penetration 
depth by far-infrared spectroscopy^^ as well as using c- 
axis YBa2Cu307/Pb Josephson junctions with a mag- 
netic field oriented parallel to the a- or 6-axisi^ Other 
studies, neutron scattering on flux- line lattices^ and 
single crystal torque-measurements )^ indicate a smaller 
anisotropy of 1.2. Related, surface impedanceS. and re- 
sistivity measurements^" demonstrate an anisotropy of 
Rs\\a/Rs\\b ~ 1-5 to 1.6 and Pal Ph ~ 1-5 respectively. 

Also, implications for the anisotropy of the supercon- 
ducting gap have been discussed. Raman scatteringli ev- 
idences a real isotropic s-wave admixture of 5%; thermal 
conductivity measurements in a rotating magnetic fieldi^ 
place a maximum of 10% based on the node positions. 
Angle-resolved photoemission spectroscopy (ARPES)-^ 
indicates larger a6-anisotropy of A^/A^ = 1.5. The 
use of untwinned single-crystals is considered crucial in 
all these studies. However, clear consensus on subdom- 
inant order parameter symmetries is not reached, nor 
is detailed angle-resolved data in the a6-plane of thin 
films available, although first attempts on twinned films 
have been performedJ^ In view of this, we present here 
new results on the anisotropy, comparing untwinned and 
twinned YBa2Cu307. 




FIG. 1: (Color online) Angle-resolved electron tunneling with 
YBa2Cu307/Au/Nb ramp-type junctions oriented every 5° 
over 360°. The YBa2Cu307 base-electrode (red) is covered 
by SrTi03 insulator (green), and contacted by a An barrier 
(not visible) and a Nb counter-electrode (blue). The arrows 
(white) indicate the main crystal orientations in the a&-plane 
of the high-Tc superconducting material. 



In untwinned YBa2Cu307 thin films, the usual 'ran- 
dom' exchange of the a- and 6-axis is eliminated. This 
enables to study the electronic properties angle-resolved 
in the a6-plane. The experimental layout is summarized 
in Fig. n Basically, the YBa2Cu307 base-electrode is 
patterned into a nearly circular polygon, changing the 
orientation from side to side by 5 degrees. A Au barrier 
and Nb counter-electrode contact each side. In this way, 
the angle with respect to the (0 1 0)-orientation is varied 
as a single parameter. 

First, bilayers of 170 nm YBa2Cu307 and 100 nm 
SrTi03 are grown by pulsed-laser deposition (PLD) on 
single-crystal SrTi03 substrates. The YBa2Cu307 films 
are optimally doped, with Tcfl > 89 K. Ramps arc ion- 
milled in the bilayers using a photo resist stencil. To as- 
sure equivalent ramp quality over 360° , the sample stage 
is rotated around the substrate normal, while maintain- 
ing the angle of incidence of the Ar-ion beam constant at 
40° with the substrate plane. The resulting ramp-angle 
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twinned 



untwinned 



SrTiOa 
YBaaCuaOr 



a//3 
a 



0.127119.0° 
0.207 97.3° 
3.866(3) 



Cell Par. (A) h 3.867(4) 



1.07°/ 357.9° 
0.75°/ 346.3° 
3.849(6) 
3.892(7) 



c 11.678(3) 11.703(7) 



TABLE I: Vicinal angle a and its orientation 13 of the SrTiOs 
substrates and YBa2Cu307 films. The YBa2Cu307 data is 
obtained after device completion. 



with the substrate plane i& = 30°. On a microscopic 
scale the interfaces may present some faceting, albeit less 
than in e.g. grain boundaries. This faceting is not ex- 
pected to affect the main conclusions of the presented 
studies. After removal of the photo-resist stencil and a 
short 90°-incidence ion-mill for cleaning purposes, a 5 nm 
YBa2Cu307 interlayer— is deposited to prepare an in- 
situ interface to a 30 nm Au barrier formed also by PLD. 
Then, a 160 nm thick Nb counter-electrode is dc-sputter 
deposited through a lift-off stencil. Special care is taken 
to obtain a clean Au/Nb interface by a 50 s rf-plasma 
etch just before Nb deposition. After lift-off, redundant 
Au and YBa2Cu307 interlayer material is removed by 
Ar-ion milling. The junctions are 4 /im wide. 

The twin behavior of (00 l)-YBa2Cu307 films is influ- 
enced by the substrate vicinal angle a and its in-plane 
orientation (3m^ Here, a is defined between the crystal- 
lographic and optical substrate-normal, and f3 describes 
the in-plane orientation with respect to the SrTiOs (100) 
crystal axis. The degree of twinning can be controlled 
from completely untwinned to the presence of four ab- 
orientations, varying a from ^ 1.1° to a small vicinal 
angle ('- 0.1°), where P = 0°. For a = 1.1°, growth 
with the 6-axis along the step-ledges is induced, and only 
one crystal orientation is present. On the contrary, 4 
twin orientations are present for small vicinal angle sub- 
strates. The twin orientations have pair-wise their in- 
plane diagonal of the YBa2Cu307 crystal aligned with 
each substrate diagonal, so that a- and 6-axes and vice 
versa are arranged nearly in parallel. After completion 
of the device, the YBa2Cu307 base-electrode is exam- 
ined with X-ray diffraction (XRD). An average of the a 
and b unit cell dimensions is found for twinned films (see 
Table P). For untwinned films, the individual a and b 
unit cell parameters can be distinguished and are close 
to single-crystal values. Detailed /i/c-scans of the (0 3 4) 
reflections show 4 different orientations for YBa2Cu307 
films grown on small vicinal angle substrates (Fig. 
associated with the above-mentioned 4 twin orientations. 
For films grown on substrates with a = 1.1° however, 
only one orientation is present (Fig.[2l3). 

The XRD and electrical data presented in this article 
correspond to the same samples. Figure 3 presents the 
electrical characterization of the twinned base-electrode 
sample (a-c), and the untwinned one (d-f,h). During the 
measurement, the magnetically shielded sample space re- 
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FIG. 2: (Color online) Logarithmic contour plots of /ife-scans 
near the (034) reflection of YBa2Cu307: (a) grown on a small 
vicinal angle SrTiOs substrate (a « 0.12°), and (b) grown on 
an Q = 1.07° and /3 = -2.1° vicinal SrTiOs substrate. Both 
scans are measured after device completion. 



duces background fields below 0.1 /iT. Trapped flux in 
or near the junctions is excluded by systematic Ic{B) 
measurements, assuring a correctly determined critical 
current density (Jc). The superconducting properties of 
the Au/Nb bilayer are independent of the orientation. 
Therefore, Jc depends on the in-plane orientation with 
respect to the 6-axis of the YBa2Cu307 crystal only, 
and presents four maxima for both samples, approach- 
ing zero in between. This is in agreement with predomi- 
nant (i3;2_j,2-wave symmetry of the superconducting wave 
function in one electrode only, and a cos(26')-dependence 
is expectediii In closer detail, the nodes of the untwinned 
YBa2Cu307 sample are found at 5° from the diagonal be- 
tween the a- and 6-axis. This presents direct evidence for 
a signiflcant real isotropic s-wave admixture. A first es- 
timate for the s- over dr^2_y2-waYe gap-ratio is calculated 
as I cos(26'o)| = 17% for a node angle 6*0 = 50°. For the 
twinned base-electrode, the nodes arc found at the diago- 
nal, which is expected if all twin orientations are equally 
present, and contributions of subdominant components 
average out to zero. 

The suppressed Jc in the nodal direction 
(< 0.01 Jc||(oio)) suggests small, if not absent, imagi- 
nary admixturesji^ for instance of isotropic is-wave or 
jrf^y-wave type which in contrast would lift the nodes. A 
significant real dxy-^s^e admixture is excluded, because 
this would induce a rotation in the same direction with 
respect to the crystal of all nodes. 

In the imtwinned case, the Jc-value is 1.8 times larger 
in the b- than in the a-direction. Preparation effects can 
be eliminated, since circular symmetry with respect to 
the substrate normal has been conserved at all phases 
of the fabrication. The normal-state resistance [RnA) is 
lower along the b- than along the a-axis, and presents a 
two-fold symmetry axis for the untwinned case. Using 
the angle-resolved values, the anisotropy in the IcRn- 
product amounts to IcRn\\b/ IcRn\\a — 1-22. 

To estimate the /ci?n-products in our junctions, we 
model them as SINS' structures, where S is YBa2Cu307, 
I is the YBa2Cu307/Au interface barrier, which has 
much higher resistance than the Au (N) and the Au/Nb 
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FIG. 3: (Color online) Jc, RnA, and Jc i?n-product vs the junction orientation with respect to the YBa2Cu307 crystal for 
twinned (a-c), and untwinned (d-f,h) YBa2Cu307 at T = 4.2 K and in zero magnetic field, (see text for description fits) (g) 
Fit parameter evolution for untwinned case, and (h) corresponding fits. 



(N/S') interface. From independent resistance mea- 
surements on our PLD Au {pAu — 4.6 pQ, cm at 
4.2 K), the mean-free-path is Iau — 18 nm, and the 
dirty-hmit coherence length is ^au — 49 nm. Using 
these values and the YBa2Cu307/Au interface resis- 
tance RbAy/Au > 10^^ Q. CTo? , the transparency at 
this interface is low: JBy/au — 440, where — 
RbA/ pAui.Au'^ From a small Fermi- velocity mismatch, 
we estimate the Au/Nb interface transparency much 
larger, JBa^/^^ < 20. The electrode-separation is cIau = 
26 nm for 30 nm thick Au and ramp angle an = 30°. 
Since Iau < dAu and Iau < S,Au, Au is in the diffusive 
regime, while YBa2Cu307 is in the clean limit with the 
anisotropic gap function Ay. 

We extend the expression for the supercurrent in diffu- 
sive SINS' structuresi^ to our case of a low-transparent 
junction between a clean d-wave superconductor and a 
diffusive NS' bilayer. The contribution of midgap An- 
dreev bound states is small in such a junction^ and can 
be neglected. 



IsRn ^ 1 1 ^'/'^X sin(x) D 7 sin(Av3) 



2TTkBT ■s—^ 

7= ^—2^ 



irkBTcNb^Nb 



T^kBTcNb + 1Bau/n, (dAu/Uu) V'^Nb + 



(1) 

(2) 
(3) 



Here x is the angle with the interface normal, and 
N = jy d(/)(ix sin(x)i? is the normalization constant. The 
integration is performed over angles (p = to 27r, and 
X = to ^ of a half-sphere of all trajectories: for each 
junction orientation, and taking the crystal orientation 
and ramp- angle into account. The barrier transmission 
coefficient D ~ cos(x) exp[K{l — cos~^(x)}] is in the 
limit of a small YBa2Cu307 Fermi velocity, where k de- 
scribes the tunnel-cone size. The sum in Eq.(2,3) is taken 
over the Matsubara frequencies u)n = 7rT(2n -I- 1). Ay 
is the anisotropic gap function in YBa2Cu3 07, and $ is 
the isotropic proximity-induced gap function in Au. A^rft 
and TcNb are the bulk pair potential and the critical tem- 
perature of Nb, respectively. The critical current Ic and 
/c-Rn-product should be found by calculating a maximum 
of IsRn over the phase difference Aiy9 across the junction. 

Tunneling along the a- and 6-axis may then be 
compared theoretically in terms of the ratio F — 
IcRn\\bl IcRn\\a- Using Eqs.(l-3), it can be shown that for 
constant properties of the Nb/Au bilayer, the ratio of the 
YBa2Cu307 gap for these directions is Ay||;,/Ay||a > F. 
Therefore, the observed anisotropy of F = 1.22 repre- 
sents a lower limit for this gap ratio, which is valid for ex- 
tremely small ratios Ay / Ajvb ^0.1. For increasing ratio 
Ay/Ajvb, the value F = 1.22 requires a rapid increase 
of Ay||ft/Ay||a. In this estimate, F depends only on the 
gap ratios and on the Au/Nb interface transparency. 

The anisotropic gap in YBa2Cu307 depends on the in- 
plane angle 9 (0 to 27r), and the angle (-^ to ^) with the 
ab-plane. Various possible symmetry functions exist to 
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describe this gap. Here, we consider the following 3D gap 
function in YBa2Cu307 consisting of a dominant d^2_y2- 
wave component with an isotropic and an anisotropic s- 
wave admixture: 

2 

Ay = Ay„ cos'i'n) c^{ cos^{0) " siu" {9)]' (4) 

with the coefficients ci > cq,C2 and ci + co + C2 = 1. 
Here Ay^ denotes the magnitude of the YBa2Cu3 07 
gap at the interface. Consistent with our earlier es- 
timate Ayj|f,/Ay||a > 1.22, the gap-ratio is taken as 
Ay||ft/Ay|[Q = 1.5 in agreement both with the observed 
node positions and ARPES.^"^ With this, the coefficients 
are found from the fit to the data Cq = 0.15, ci = 0.83 
and C2 = 0.02. Other choices for the gap symmetry func- 
tions lead to slightly different numbers but will not alter 
the basic results of our calculations. 

A series of fits is presented in Fig. Oi: the wider the 
tunneling cone (smaller k), the smaller is the width of 
the oscillations in the /c-Rn(6')-dependence (arrow). The 
effective YBa2Cu307 gap Ay^ and 7b^„/„6 must then 
become larger. This dependence is presented in Fig. 
The minimum value Ay^ occurs for normal-incidence tun- 
neling (limK — > oo), so that Ay^ > 6.4 meV. For reason- 
able Ay^ values (< 0.5 eV), IBA^/Nh varies from 13.6 to 
about 18. This gives an estimate for the Au/Nb interface 
resistance: RbA = 0.36 ± 0.05 nH. cm^. In contrast to 
"fBAu/Nb^ it is not possible to give an accurate estimate for 
the YBa2Cu307 gap from our data. Therefore, we choose 
to fit the data with Ay^ = 44 meV, 7b^„/„6 = 16.8 in 
Figs.|3t and|3f. These are not claimed to be the correct 
values; the simulation demonstrates however that large 
Ayj, may well be consistent with small IcRn values. For 
the untwinned case k = 18.4, corresponding to a tunnel 
cone with a full- width-half-maximum (FWHM) of 31.0° 
(cosine-term of D not included). For the twinned case, 
the /ci?n (0)-dependence is simulated with the same pa- 
rameters, except for a slightly smaller cone (k — 26.3, 
FWHM = 26.0°), and assuming equal presence of both 
twin orientations: i[/ci?„(0) -l-/ci?n(^+ f )]■ The smaller 
tunnel cone for the twinned case is consistent with higher 
RnA- and lower Jc-values. This may result from a 
slightly thicker tunnel barrier at the YBa2Cu3 07/Au in- 
terface, e.g., due to minor variations in the Au PLD- 
conditions, modifying this interface. 

The i?,„A((?)-dcpendcncc is fitted with an ellipsoidal 
relation of the conductivity projections along the main 



crystal directions of the YBa2Cu307. Written in terms 
of the i?„A-values along these directions, this gives 

RnAi9) = y^i?„Ajj, cos2(0) -|- i?„A|^ sin" (6). Fig. H 

shows the result using i?„A||Q = 44 nft cm^ and 
= 29 nQ cm^. For the twinned case, a 
geometrical average of the conductivities is assumed, 
RnA{9) = 2/[Rr,A-\e) + RnA-^9+^)]. The used val- 
ues in Fig. read i?„^||a — 141 nQ cm^, = 
47 nil cm^ . Although these phenomenological fits are in- 
dicative, angle-resolved calculations including aspects of 
the YBa2Cu307 band-structure and band-bending effects 
are needed for a detailed understanding. Finally, the 
Jc{9) fits are obtained with the ratios of the IcRn{9)- and 
the _R„^(0)-dependencies, the ensemble of which gives a 
consistent simulation of the angle-resolved junction prop- 
erties. 

The experimental results support theories based on 
a 2-band model of the chains and planesSLS^ modeled 
with a symmetric, anti-symmetric and isotropic compo- 
nent. Furthermore, our findings agree with c-axis tun- 
neling from two twinned YBa2Cu307 grains to a Pb 
counter-electrode that depends on the magnetic field 
orientation,^" and angle-dependence studies on grain 
boundary junctions.^* For all-high- Tc junctions and cir- 
cuits, we mark the anisotropy as a possible intrinsic 
source of their limited reproducibility: both twin orien- 
tations may not be uniformly present, yielding an impor- 
tant variation in J^. Control over the crystal orientation 
then presents a key to improvement. Another important 
aspect concerns the nodes at 5° from the (1 1 0) crystal 
direction. The best choice for the electrode-orientation 
of devices aiming a cZ-wave induced second harmonic in 
the current-phase relation, such as ^-SQUIDs based on 
grain boundary junctions, may therefore deviate from the 
(110) crystal direction. 

In conclusion, an angle-resolved electron tunneling 
study using Josephson junctions with an untwinned 
YBa2Cu307 base-electrode is presented. Evidence for 
significant in-plane anisotropy in the electronic proper- 
ties of YBa2Cu307 is found. 
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